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1. Introduction 



For the description of the vacuum creation of an electron-positron plasma (EPP) under the 
action of a strong time dependent external electric field (the dynamical Schwinger effect) the reli- 
able tool is the kinetic equation (KE) It is an exact non-perturbative consequence of the basic 
equations of motion of QED for a linearly polarized spatially homogeneous electric field ("laser 
field"). In the present work we continue to analyse a quadrature solution of this KE obtained in the 
low density approximation |^. We estimate here the role of photon cluster processes introduced 
1^] in comparison with the orthodox multiphoton mechanism and show that this mechanism plays 
the dominant role in the case of a large photon number absorbed from the external field reservoir. 
It corresponds to the optical and X-ray range of the laser radiation. The quasiparticle EPP excita- 
tion within a period of the laser pulse action was investigated rather detailed (e.g., [||]). Below we 
discuss the mechanisms of the residual EPP generation. 

We restrict ourselves here to the consideration of the multiphoton domain corresponding to 
large adiabaticity parameters 7= {E^-/Eq) ■ (v/m) ^ 1 and subcritical fields Eo<^Ec = m^/e . The 
basic equations are represented in Sect. 2. Corresponding to the multiphoton mechanism of the 
EPP excitation the nonlinear harmonic analysis is discussed in Sect. 3. Sect. 4 is devoted to an 
investigation of the momentum spectrum of the problem and the threshold conditions of the energy 
gap breakdown. The proof of the dominating role of the photon cluster process in the domain of 
large total photon numbers is given in Sect. 5. The distribution function in this approximation is 
obtained here also. The basic results of the work are summarized in Sect. 6. 



2. Low density limit 

For the description of the EPP vacuum creation we use the exact nonperturbative KE obtained 
in the work for an arbitrary time dependent and spatially homogeneous electric field with linear 
polarization. In Hamiltonian gauge, A^{t) = (0,0,0,A(?)), the electric field strength is E{t) = 
—A{t). The one-body electron (positron) phase space distribution is 

t 

/(p,0 = ^A(p,0 1 Jf'A(p,O[l-2/(p,O]cos0(?,O, (2-1) 

to 

where 

X{p,t)=eE{t)eJ(o\p,t) (2.2) 



denotes the amplitude of the EPP excitation, co(p,f) = y£^(p) + (/7|| — eA{t)y stands for the 
quasienergy with the transverse energy e± = (m^ + p\y^^, and the high frequency phase is 

e{t,t')=2 j'dxco{p,x). (2.3) 

The distribution function in the quasiparticle representation is defined for the in-vacuum state, 
/(p,f) = (/?i|a+(p,f)a(p,f)|/«). For the generalization of the KE (2.1) to arbitrary electric field 
polarization see [|5|]-[^. 
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In the low-density limit, / ^ 1, KE ( |2.1[ ) obeys the formal solution 

2 

1 



/low(P,0 



(2.4) 



were the external field is supposed to be switched on in the infinite past, — )• — oo. The low- 
density approximation implies that the external field is weak, i.e., E <C Ec- Below we will restrict 



ourselves to the analysis of Eq. (2.4) in the limit f — oo, which defines the momentum distribution 
of the residual EPP [3], 

1 ' 



foutiv) = lim//o„.(p,0 



dtX {v,t)e 



(2.5) 



(To arrive at ( p3| ) the representation of the phase ( P3| ) via antiderivatives has been used, d{t,t') = 

e{t)-d{t').) 

Hereafter we will consider the domain of the multiphoton mechanism of EPP vacuum creation, 
where the adiabaticity parameter defined above is large [8], 7^ 1. A perturbation theory in the 
small parameter 1 // ^ 1 can be constructed here. In the case of a monochromatic laser field 



= (£'o/v)cosvf, £■(?) =£'osinv?, 



this statement is based on the representation of the quasimomentum 

m 



P = P\\ 

As a first step let us consider (|2.5[) in an expansion. In leading order, one obtains 



— cosvf. 
7 



(2.6) 



(2.7) 



A(p,?) -> A()(p)£(0, A()(p) =e£x/w^(p), 



(2.8) 



and d{t) = 2(0ot. This means that contributions of the high frequency harmonics of the external 



field are neglected in Eq. (2.5), yielding the result 



fout{p) = \^\E{(0 = 2(Oo)\\ 



(2.9) 



where E{(o) is the Fourier transform of the field strength E{t). In this approximation, the vacuum 
EPP production takes place when for the frequency of the one-photon e^e^ pair creation process 
(Oiy holds (Oiy = IcOq. This mechanism is exclusive here. Its intensity is regularized by the presence 
of the frequency coiy in the spectrum of an external field. 



Let us consider the special case of a monochromatic field (2.6). The distribution function /(p) 
becomes 



1 



foutip) = -Xe'o8 [2(Oo{v) - V] 5(0) 



(2.10) 



Substituting 5(0) — ^ 8{lK/Tp) — ^ Tp/ln (where Tp ^ 2k/v is the period of the field pulse action), 
we obtain the spectral density for the production of the EPP 



dt 



1 



— Ao%^5[2o^,(p)-v] 



(2.11) 
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The 5 distribution defines iiere the admissible momentum set for v ^ 2m. Thus, employment the 
field ( ^ ) nonconvergent for t — > ±00 leads to the singular distributions ( 2.10 ) and ( 2.11 ), for which 
the physical meaning is restored by calculating physical quantities such as densities. 

In order to facilitate breakdown of the energy gap by EPP creation, one can take into ac- 
count the multiphoton mechanism of EPP production. These processes are described by the high 



frequency multiplier in Eq. (2.5). The relevant methods will be developed in the next Section. 



3. Multiphoton processes 

Let us use the non-perturbative method of photon counting assuming that the electric field 
is periodical A{t) = A{t + 2n /v), where V = In/T is the angular frequency and A{t) = A{—t), 
in order to provide the property of being an even function of the quasiparticle energy, ft)(p,f) = 
ft)(p, — f). The corresponding Fourier transform 



^(P)0 = ^ n„cosnvf 

n=0 



leads to the decomposition of the phase in Eq. ( p.5| ) 

B{t) = 2Q.Qt + 2^ a„sin/ivr 



(3.1) 



(3.2) 



where a„ = 2n„ /vn . In case of a harmonic external field, Q.q is the renormalized frequency [p 



Let us employ now the approximation ( |2.8| ) and use the non-perturbative decomposition 

00 

exp(/asin0)= £ Jkia)e""^ , (3.3) 
where Jk{a) are the Bessel functions of order k. Let us now consider the integral in Eq. (p3|) and 



perform the substitutions (3.1)-(3.3), leading to 



f dtE{t)e^''''>' j/o(2no) + n I J-k{an)e-'^'"'^'" + ^ -^)| . (3.4) 



The first term with /o(2^^o) describes the direct vacuum excitation at the frequency (Oiy = 
2Q.0 and corresponds to Eq. ( [2. 10| ). If we are interested in excitations at lower frequencies, this 
contribution can be omitted. One can omit in Eq. (3.4) also the series with the substitution k ^ —k 



since the corresponding frequencies cannot lead to a compensation of the high frequency phase 
2^0^- At last, let us replace the infinite series and products by finite ones to arrive at the expression 



Npi, 



Jip,Nph,Nc) = / dtE{t)e^'^' n L J-kMe 



-i(kn)vt 



(3.5) 



«=u=i 



Here the multiphoton processes of the energy absorption from the photon reservoir of the 
external field are characterized by the pair of indices n and k. The index n corresponds to the 
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ordinary multiphoton process (see its source in Eqs. (3.1), (3.2)), whereas the index k marks the 
group (cluster) multiphoton process when each of k clusters contains n identical photons with the 
energy v. The ordinary multiphoton process corresponds to the simplest "cluster" of the order k=\. 



The appearance of photon clusters in the multiphoton processes in Eq. ( |3.5| ) is a consequence of the 
nonlinear field dependence of the quasiparticle energy G)(p,?). The probability for the generation 
of a n-photon cluster is defined by the amplitude a„ in the decomposition (|3.2[) (the argument of 



the Bessel function 7/t(a„) in Eq. ( |3.5[ )) while the probability of a simultaneous /c-cluster absorption 
corresponds to the Bessel function of the order k. 



Now one can identify the numbers Npi, and Nc in Eq. (3.5): Nph is the maximum number of 
photons in a cluster, and Nc is the maximum number of clusters. The limit Nph °°, Nc ^ °° 



corresponds to the exact formula (3.4) 



In order to perform the integration over t, it is necessary to carry out the summation in Eq. (|3.5D. 



For this aim let us remark that products of the different series in Eq. (3.5) among themselves 



generate the exponent product. Thus, for the simplest harmonic field (2.6) we obtain 



J{p,Nph,Nc) = -imt.{^i^^o-vk^ph + v)- 5(2^0 -vk^pH-v)}llJ-kian) ,(3.6) 

k=l n=0 



where the total photon number in a cluster of first order (k = 1) is 

Npi, ^ 

^ph=Y."= 2^ph{Nph + 1) . 
n=l 



(3.7) 



Let us substitute now the relation (3.6) into Eq. (2.5). This leads to a sum of products of two 
5-distributions. They refer to (i) the energy conservation law in the multiphoton process 



2ao(p)-vA^,^, = 0, 



(3.8) 



where Nj^i = kJ/'ph ± 1 is the total number of the photons taking part in the process (contribution 
±1 is stipulated by the multiplier E{t) in Eq. (^^)), and (ii) to the argument of the cluster processes 
in the product J J* leading to either k\=k2 = k ox {k\ — k2)^ph = i2. The last condition can not 



be fulfilled for any integer numbers k\ and ^2- Thus, we obtain for the distribution function ( [2.5| ) 

foutip,Nph,Nc) = 
Here, we used the relation 



— yi£{5[2ao-vA^+] + +5[2ao-vA^,,]}n/-.M • (3-9) 



n=0 



Y,5[v^pH{k-k')]Fu^y 



1 



V^ph 



■Fk.k 



(3.10) 



for an arbitrary function F]^^]^!. Such an approach allows to avoid the hypothesis of a phase random- 
ization [3]. 



4. Spectrum 

Our next goal is the investigation of the equation (p.8[), which defines the threshold photon 



numbers and the spectrum of the created EPP. For a given frequency v, the conditions (3.8) fix the 
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minimum combinations 



2ao(p = 0)/v±l 



(4.1) 



of the threshold photon {N'^i) and cluster {N'J) numbers, whereby the breakdown of the energy 
gap is possible for the produced EPP. Absorption of a greater photon number with k^ph > N'/N'^j^ 
generates an EPP with p ^0. 

The first step in the analysis of the threshold conditions is the determination of the momentum 
dependence of the renormalized quasienergy no(p) in explicit form. According to Eq. (3.1 ) 

1/2 



^O(P) = ^ 



7C/V I 

(o(n,t)dt = — 

-jz/v 71 



du ■ 



el + {p\\-um/YY 

1-m2 



(4.2) 



for the field (|2.6|). From here it is possible to receive various approximations. The upper bound is 

no(p) = 



com 



■2p\\m/y , 



(4.3) 



where (p) is the renormalized quasienergy with the mass = m^J\ + \/y^ being renormalized 
in the field. From Eq. (43) one sees that the anisotropy effect is proportional to I/7. In the isotropic 
approximation we obtain Ilo(p) = w*(p). 

From Eqs. (3.8) and ([4.3[) we find the momentum spectrum of the created EPP 



Pl,2/'«* = ^COS(/) lb 



n 2 



1 + 4rCOS^0 
f 



1/2 



(4.4) 



where is the angle between the vector p and the axis p^ = p\\. 

Let us consider below the isotropic approximation. The momentum spectrum ( |4.4j ) transforms 
then according to 



1/2 



1 



(4.5) 



For p = Q the threshold condition follows from here (the symbol [x] means the integral part of x) 

1 



^^N';N%{N'J + l)±\ = [2m,/v], 



(4.6) 



which defines the set of the threshold numbers N'J" and N'^f^ for the given total photon number 
that is necessary for a breakdown of the energy gap. For N^g, = 1 , the one-photon mechanism of 
EPP excitation is realized again. Combinations of the numbers N'J and A^^^ (compatible with the 
threshold condition (4^6)) are presented in Table 1 for the simplest few -photon processes with the 
minimum number A,^,. From this table one can see that number of the photon degrees of freedom 
is larger than for conventional multiphoton processes, at least in the domain of the few-photon 
mechanism, i.e. for v ~ m. 

Let us introduce now the superthreshold photon {5n) and cluster {5k) numbers. The corre- 
sponding total numbers will be equal Nph = A^'/, + 8n and Ac = N'J + 5k. The numbers 8n and 8k 
define the discrete momentum spectrum 



PSnM = 



Am, 



-5{N„Nph) + \ 



1/2 



(4.7) 
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f(0.P3) 




Table 1: Distribution of the total photon 
numbers A^tot to clusters 
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Figure 1: A typical distribution function of the resid- 
ual EPP in periodic field with Gaussian shape (E = 
Q.2Ec,X ^Q.Qlnm) 

where 

5iNc,Npi,) = N'Jd^ph + ^;i5k + 8^ph5k . (4.8) 

Thus, we obtain a discrete spherical symmetrical (in the isotropical approximation) momentum 
spectrum. For v <^m and 8n <^ N'^j^ and 5k <^ N'J, one obtains 

PSnM = {^*v^{N,,Nph)Y'^ , (4.9) 

i.e., the spectrum is quasicontinuous in this case. 

Using the solution of Eq. ( |3.8[ ), we can rewrite the 8 distributions in Eq. (3^) in the form 



5 [2Q.0 - vN, 



±1 

tot\ 



-5ip-ps„,Sk) 



(4.10) 



^PSn,Sk 

The next subsidiary problem is an analytical estimate of the Fourier coefficients of the series 



(3.1) 



jz/v 



n > 1 



(4.11) 



dt ft)(p,?)cos«v?, 

71 J-nlv 

and the corresponding arguments a„ of the Bessel functions. It can be shown that each Fourier 
coefficient ( |4.1l[ ) can be represented by a series expansion with respect to the adiabaticity parameter 
1/7 ^ 1- The minimum power of such a kind of decompositions is fixed by the harmonic number, 
~ l/y"- It is important that the even harmonics only correspond to the isotropical distribution 
in the momentum space. As a result we obtain for the even harmonics 

n—\ ( m 
~ n-«!22«-2 

Using the definition A^^f, Eq. (3.7) and Table 1 (for small numbers n and K), we obtain for Ho = ft) 



V 



25', s= 1,2, ... 



(4.12) 



V 



kn{n+l) ± 1 



(4.13) 



and hence af^l for large photon numbers and a,f ^ 1 for small ones. This allows to use the 
asymptotic representation of the Bessel functions for small or large values of the arguments. 
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5. Role of the photon cluster processes 



The following working formula follows from Eqs. (3.9) and (4.1C) 



fout{p;Nph,Nc 



^{p-psnM) n -^h 

11=0 



k=0 



(5.1) 



Here, the even harmonics are considered only. 



In order to use the formula (5.1), it is necessary to find the threshold numbers N'J, and N'/. 



Equation ( |4.6[ ) serves for the definition of them. For a fixation of these numbers the following 
prescription is offered: for a given frequency v the different orders N'/ of clusters are tested and 
for each number A'^^'.'" the corresponding photon number N^j^ is found. These threshold numbers 
correspond to the creation of an EPP where all pairs are at rest. The superthieshold numbers with 
5n^0 and 5k ^0 describe the moving EPP. 

Now we compare the efficiency of the usual multiphoton and cluster mechanisms. The sim- 
plest method is a comparison of two limiting values of the distribution function for = corre- 
sponding to either the maximum photon number in one cluster (the first limiting case, the ordinary 
multiphoton process) or the maximal cluster number with the minimal photon number (n = 2) in 
each cluster (the second limiting case) with the same total threshold photon number in both cases. 
Let us assume that the total photon number is large 



ph 



Then we have (i) in the multiphoton case 



1, Nu 



and (ii) in the multicluster case 



N'pf, = 2 (the second harmonic only), N'J = Ntot/'i = [2m*/3v] . 



(5.2) 



(5.3) 



(5.4) 



Basing on Eq. (^), we rewrite the distribution function for the threshold values of the photon 
numbers as 



■M'-i^N';,)S{p) . 



(5.5) 



2v.yVph(0^ p 

In order to avoid the infrared divergence problem, we introduce the clustering coefficient de- 
fined as the ratio of two limited distribution functions: the multicluster (/c) and multiphoton (fph) 



distribution functions with the corresponding threshold numbers ( p.3| ) and ( |5.4[ ) 

^ = fc/fph ■ 



Substitution of Eqs. (4.12) and ( p3| ) leads to the strong inequality 

^ >1 . 



(5.6) 



(5.7) 



Thus, the photon cluster mechanism is dominating in the region of large total photon numbers. 
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This conclusion allows to generalize Eq. ( |5.5| ) for the case of an EPP with pairs of nonzero 
momentum, that corresponds to the superthreshold photon cluster numbers Nc ^ N'J , 



■PSn,Sk) 



k=0 



where according to Eq. ( ft.l2[ ) 



mz 



a2 



(5.8) 



(5.9) 



The structure of Eq. (5.8) ensures a slow power (nonexponential) decrease of the distribution func- 
tion in the momentum space. This distribution is very extensive (of the order of m) and rapidly 
oscillating relative to some average value. 

The transition to the superthreshold cluster numbers Nc » N'/ in Eq. (5^) is reflected here in 
two effects: the appearance of the higher cluster harmonics (index of the Bessel functions) and the 
account of EPP motion (p ^ 0). The first feature describes, apparently, the foliated structure of the 
distribution function [||] whereas the second one defines the momentum dependence on it. Below 
we will consider the last effect only. Then, in this approximation, we obtain from Eq. (5^) in the 
limit Nr ^ °° 



foAp) = '?^t4M)5{p-p') , 



6v«3 p 



(5.10) 



where according to Eq. ( |4.10| ) p' = 2k that corresponds to two photons in a cluster. Substitution of 
the summation by an integration leads to the result 



(5.11) 



According to Eqs. ( fl.l3| ) and ( |5.9| ), the relation I holds. Using the corresponding asymptotic 
representation of the Bessel functions 



Jniz) = y^2/7lz COS (z - Kn/2 - k/A) , z»l, (5.12) 
one can obtain from Eq. ( |5.1l| ) the distribution function in the considered approximation as 

fou,{p)=fZ{p)+fZ{p), (5.13) 



where 



is the smooth part and 



'icolpml 



/r.(p) 



fZi9) = fZ{v)-<^o, 



mz 



872 w*V 



■Tl{n+\/2) 



(5.14) 



(5.15) 



is the oscillating part in the momentum space. 

Thus, the simplest approximation (5.13)-( ^?T^ ) reproduces reasonably well some important 
features of the results of numerical calculations of the distribution of the residual EPP on the basis 



of Eq. (2.5), see Fig. 1. 
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6. Summary 

In the present work, we have developed further the ideas of ^ which initiated the consid- 
eration of the dynamical Schwinger process as a multiphoton process. Here, the new class of the 
multiphoton processes to wit the photon cluster one was introduced. Such a mechanism describes 
an other analytical tool in comparison to the orthodox multiphoton processes and is stipulated by 
the nonlinear field dependence of the quasienergy. The photon cluster process is a new class of 
cooperative effects in which the photon set behaves as one photon with the total energy of the set. 

We analyzed the distribution function of the residual EPP in the multiphoton region in the low 
density approximation. The basic result of our work is proving the statement that the photon cluster 
process is dominating if the total photon number is large. Our preliminary analysis has shown also 
that the orthodox multiphoton process plays the main role in contrast to the case of small total 
photon numbers. 

We had obtained the final formula for the distribution function of the residual EPP in the region 
of action of the cluster mechanism EPP excitation and noted some its features: 

- the distribution function has cylindrical symmetry and is flattened in the direction of the 
acting electric force; 

- the distribution function is very extensive (of order m) and rapidly oscillating in the momen- 
tum space. 

Subsequent development of this research direction addresses has not only the current prob- 
lems (calculation of the number density EPP, problem of the infrared divergence, estimation of 
the distribution function in the dominant domain of the orthodox multiphoton mechanism etc), 
but also more deep questions as well (e.g., understanding the character of the rapid oscillations in 
momentum space (Fig. 1), correlation between the quasiparticle and residual EPP). 
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